Technical  Report  JSR-74-6 


June  1975 


LOW-FREQUENCY  SOUND  PROPAGATION 
IN  A FLUCTUATING  INFINITE  OCEAN-11 


By  F.  ZACHARIASEN 


Th#  viiw  and  conclutiom  contained  in  this  document  aia  thoae  o*  the  authon  and  thould  not  be 
interpreted  ai  nacetsanly  representing  the  otfic.al  policies,  either  expressed  or  implied  o*  the 
Defense  Advanced  Research  Protects  Agency  or  the  U S.  G'  -rnment 


Approved  for  public  release;  distribution  unlimited. 


STANFORD  RESEARCH  INSTITUTE 
Manlo  Park,  California  0402S  • u.S.A. 


Copy  No. 


Contract  No.  DAHC15-73-C-0370 

ARPA  Order  No.  2504 

Program  Code  No.  3K 1 0 

Date  of  Contract  2 April  1973 

Contract  Expiration  Date:  30  June  1975 

Amount  of  Contract.  $1,297,321 


DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 
ARPA  ORDER  NO.  2504 


Sponsored  by 


UNCLASSIFIED. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whon  Data  Entorodl 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

2 GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

1 

U 

1 

w-7 -Tt  

//'.LOW-EREQUENCY  ^OUND  PROPAGATION  IN  A FLUCTUATING 

f Infinite  oc  bandit,  * i 

O'.  *P«E  OF  REPORT  & PERIOD 

COVERED 

^Technical  ^p^t/^ 

: _ 5 i 

^ri  ^■^^►3poo'|u" ' 

NUMBER 

1 

j 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Stanford  Research  Institute 
Menlo  Park,  California  94025 

io,  prograkLelemen 

AREA  '«?^trhK  UJtl.1 

T,  PROJECT,  TASK 
T NUMBERS 

l 

REPORT  DATE 

13.  NO.  OF  PAGES 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

(£puneai075j 

36 

I 

I 

| 

15.  SECURITY  CLASS,  (of  this  report) 

Unclassified 

I 

14.  MONITORING  AGENCY  NAME  & ADDRESS  (If  dlff,  from  Controlling  Office) 

| 

• 

15a,  DECLASSIFICATION /DOWNGRADING 
SCHEDULE 

/ 


ZacH i riasen 


—CONTRACT  OR  GRANT  NUMBER(s) 

w 


DAHC15-73-C-p37dh 

| 0 ¥ P 1\  \ Q \ dev-  0 1 


16.  DISTRIBUTION  STATEMENT  (of  this  report) 

Approved  for  public  release;  distribution  unlimited 


17,  DISTRIBUTION  STATEMENT  (of  tho  abstract  entered  in  Block  20.  If  different  from  report) 


IB.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  side  il  necessary  and  Identify  by  block  number) 

Internal  waves 

Sound  propagation  fluctuations 


^ABSTRACT  (Continuo  on  reverse  sldo  if  necessary  ond  identify  by  block  number) 

An  approximation  method  for  calculating  fluctuations  due  to  internal  waves  in 
sound  propagation  in  the  ocean  is  outlined,  including  the  effects  of  the  sound 
channel.  Results  are  presented  in  a simple  form  in  which  the  geometrical  optics 
limit  is  transparent, 


DDri473 


FORM 
1 JAN  73 

EDITION  OF  1 NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Doto  Entered) 


33A  £~M>- 


e: 


;/ 


\ ..  •>.  - . r.J  - * 


4 


-i 


.T'" 


I 


v, 


"i 

3' 


1 


i 


i 


H 

•i 


l 


m 


1 

II 


m 


m 


S i.H'.'ScMviS  ’ 

t ' . i t f.  iu  f(  .v‘.  -«  •.  :r.  • 


> 


CONTENTS 

DD  Form  1473 

I INTRODUCTION 1 

II  REVIEW  OF  SUPEREIKONAL  APPROXIMATION 2 

*■' 

III  FLUCTUATIONS  IN  A HOMOGENEOUS  ISOTROPIC  OCEAN  9 

IV  RYTOVS  METHOD  IN  AN  INHOMOGENEOUS  OCEAN 15 

V FLUCTUATIONS  IN  AN  INHOMOGENEOUS  OCEAN 20 


vt  rnNrT.iismM  oa 


■ 


I INTRODUCTION 


In  an  earlier  report,  we  described  an  approximation  scheme,  called 
the  supereikonal  approximation,"  for  calculating  fluctuations  produced 
by  a specified  spectrum  of  sound-speed  fluctuations,  in  pressure  signals 
from  a CW  source  in  the  ocean.  In  that  report,  our  analysis  was  confined 


to  the  case  where  the  fluctuations  were  superimposed  on  a homogeneous 


isotropic  ocean. 


Here  we  should  like  to  extend  and  improve  on  the  earlier  version  in 


two  major  ways.  First,  we  should  like  to  rewrite  the  results  obtained 


in  Ref,  1 for  the  homogeneous  case  in  a much  simpler  form,  in  which  the 
geometrical  optics  limit  valid  for  short  ranges  is  transparent,  and  in 
which  the  corrections  to  this  limit  that  became  important  at  large  ranges 


are  easily  visualized.  We  should  also  like  to  derive  formulas  for  thi 


fluctuations  in  jhe  quantities  of  most  immediate  experimental  interest — 
namely,  the  mean  square  phase  of  the  received  pressure,  and  the  mean 
square  logarithm  of  the  amplitude  of  the  received  pressure. 


Second,  we  should  like  to  extend  our  treatment  to  the  realistic  case 


in  which  the  background  ocean  is  not  isotropic  and  homogeneous,  but  con- 
tains a sound  channel.  As  we  shall  see,  we  can  handle  this  situation  as 
well,  with  only  a modest  increase  in  complexity,  as  long  as  the  fluctua- 
tions and  ranges  are  such  that  separate  ray  paths  are  still  identifiable. 
Our  result  will  then  apply  to  fluctuations  in  the  signal  associated  with 


a given  ray  path. 


Numerical  results,  and  specific  models  of  fluctuation  spectra,  will 


be  treated  in  a separate  report. 


References  are  listed  at  the  end  of  the  report. 
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II  REVIEW  OF  SUPEREIKONAL  APPROXIMATION 

Let  us  begin  our  discussion  by  neglecting  the  effects  of  the  sound 
channel.  The  problem  of  sound  propagation  in  the  presence  of  fluctuations 
superimposed  on  a homogeneous  isotropic  background  is  easier  to  set  up  and 
to  visualize  than  when  the  fluctuations  are  superimposed  on  an  inhomogo- 
nous  background,  so  it  is  conceptually  advantageous  to  work  out  this  case 
first.  Then,  as  we  shall  see  later,  when  the  inhomogeneous  background 
lepiesenting  the  sound  channel  is  introduced,  the  analysis  can  bo  carried 
out  very  much  as  in  the  homogeneous  case,  and  the  resulting  formulas, 
while  geometrically  more  complex,  are  entirely  analogous  to  those  obtained 
in  the  simpler  example. 

Our  analysis  will  be  based  on  the  supereikonal  approximation,  and  it 
will  be  convenient  at  this  point  to  give  a brief  review  of  the  results  of 
this  method.  Ihe  problem  is  to  evaluate  the  pressure  p(x)  at  a point  x 
produced  by  a point  source  (which,  for  convenience,  we  can  take  to  have 

unit  strength)  located  at  the  origin.  The  sound  propagation  from  the 

— ► 

souice  to  the  point  x is  through  an  isotropic  homogeneous  ocean,  in  which 
the  sound  speed  is  c,  on  which  is  superimposed  a fluctuation  in  sound 
speed  6c(x)  that  is,  of  course,  very  small  compared  to  c.  Mathematically, 
then,  the  pressure  satisfies  the  wave  equation 

(V  + k ) p (x)  = V(x)p(x)  (2.1) 

where  k = uj/c  for  a source  emitting  sound  of  frequency  0),  and  where 

V(x)  = 2k2  6c(x)/c  . (o  2) 
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The  boundary  condition  associated  with  Eq . (2.1)  is  that  a 


as  x - 0, 


p(x) 


4xx 


(2.3) 
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(we  shall  denote  |x|  simply  by  x) 

(We  may  comment  that  in  the  case  of  an  inhomogeneous  background, 
in  Eq.  (2.1),  k is  simply  replaced  by  k(x)  = w/c(x),  where  c(x)  is  the 
sound  speed  in  the  inhomogeneous  background.) 

Equation  (2.1)  may  be  case  into  integral  form  through  the  use  of  the 
outgoing-wave  Green's  function 


2 2 3 
(V  + k ) A(x  - y)  = 6 (x  - y) 


« 


I 
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(2.4) 


explicitly,  ive  have 


Mx)  = -i-elk  * 

4 XX 


(2.5) 


Then  we  may  write,  in  place-  of  Eq.  (2.1), 


p(x)  = A(x)  + Jd3y  A (x  - y)V(y)p(y) 


(2.6) 


Iteration  of  this  integral  equation  generates  the  perturbation  series  for 

p(x),  which  is  more  convenient  to  write  in  Fourier-transformed  form,  as 
follows : 


P (q ) = A(q)  + A (q) 


f d % 

i)  / ? V(q 

J (2x) 


)A(q  - q^ ) 


(2.7) 


+ A(q) 


r«\  r«%  _ 

J^?J^  v<q> 


)A(q  - 


qi)V(q2)A(q  - - qg) 
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whore 


A(q) 
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and,  of  course,  where 


p(q) 


/ 


3-> 

d x c 


->  -> 

-iq  .x 


p(x) 


(2.8) 


(2.9) 


The  supereikonal  approximation  now  consists  of  neglecting  all  momen- 
tum transfer  correlations  in  the  perturbation  series.  That  is,  we  approx- 
imate (q  - q - q - ...  - q )2  - k2  (-  ie  by  q2  - 2q  • (q.  + qn  + . . . + q )2 
x w 11  l z n 

+ qj_“  + c\2^  + ...  + q^2  - k2  + is,  and  neglect  all  terms  of  the  form 
q ^ * 9^  when  l t j . Note  that  the  first  approximation  occurs  in  the 
second-order  term  in  V.  Once  this  simplification  is  made,  the  perturba- 
tion series  can  be  summed  exactly,  and  one  obtains  the  result 


i [0 k2+  —■  + B I (B  ,x)+ie] 
4p 
e 


who  re 


(2.10) 


I (0  ,x) 


2 

-i [sq . x+Bs ( l-s)q  ] 
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(2.11) 


This  expression  constitutes  the  supereikonal  approximation  to  the  pressure. 


The  conditions  under  which  it  is  valid  are 


'fcliS^.TVskv-,-  «*“  ■vw&err.  toiktittw  t*K*v«r3ai£*  *<(i  • ; 


and 


k x » 1 


k L » 1 


X < 


k2L  \6c> 


(2.12) 


% 
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On  fact,  the  last  condition  may  well  be  too  stringent.)  Here  L is  the 
correlation  length  of  the  sound-speed  fluctuations-!. c.  , the  correlation 
function  C (x  - y)  = <V(x)V(y)>  vanishes  when  |x  - y|  l. 

It  is  worth  noting  that  if,  in  Eq.  (2.11),  the  0s(l  - s)q2  term  is 
omitted  from  the  exponent,  we  obtain 


p(x) 


_l  i^H1 

- e ' 0 

ditx 


V(sx)ds  lx 


Which  is  the  conventional  WKB,  or  eikonnl,  or  geometrical  optics,  approxi- 
mation to  the  Pressure.  The  primary  virtue  of  the  supereikonal  form, 
therefore , is  that  it  contains  as  limiting  cases  both  the  conventional 
eikonal  and  and  complete  first-order  perturbation-theory  approximations. 

tthile  Eqs.  (2.10)  and  (2.11)  do  constitute  a closed-form  solution  for 
the  pressure,  the  expressions  are  still  a bit  unwieldy,  and  further  simpli- 
fication of  them  is  useful.  To  this  end,  let  us  evaluate  the  integral  in 

Eq.  (2.10)  by  stationary  phase,  keeping  in  mind  that  x and  k are  both  large, 
The  stationary  phase  point  is  0 where 


k2  - t l(Box)  ♦ Bo  ± 10,3)1  - 3/2 9„  = 0 

0 0 


40 


0 


From  Eq.  (2.11),  we  may  estimate  that 


0 d3 


I(P,x) 


0=3 


6c 

c 


hence,  if 


x < 


L2  — 

6c 


(2.13) 


the  stationary  phase  point  is  accurately  given  by  the  solution  of  the 
simpler  equation 


2 

k 


2 , n 2 

x /4B0 


0 


and  is  located  at  3Q  = x/2k.  Thus  we  find 


(2.14) 


To  approximately  evaluate  the  integral  I in  this  expression,  we  recall 
that  the  supereikonal  approximation  should  be  exact  to  first  order  in  V. 
This  l gcj u irenie nt  then  yields  the  result 


- f 


P(x)  = A(x)  exp  — - J d~y  A(x-y) V(y)A(y) 
A(x) 


(2.15) 


This  expression  is  known  as  Rytov's  approximation  to  the  pressure.  A 
direct  derivation  of  it  may  be  made  by  replacing  the  wave  equation  (2.1) 
by  an  equation  for  log[p(x)/A(x)]  and  solving  this  to  first  order  in  V.2 
However  the  justification  for  the  approximation  is  somewhat  obscure  in 


6 


this  direct  derivation;  in  the  approach  via  the  supcreikonal  technique, 
what  is  being  left  out  is  more  clearly  visualized. 


In  any  event,  depending  on  the  validity  of  the  criterion  (2.13),  one 
may  use  either  the  supereikonal  expression  (2.10)  or  the  Rytov  expression 
(2.15)  to  proceed  further.  Wo  shall  use  (2.15). 

Let  us  write,  then, 


p(x)  - A(x)  e' 


(2.16) 


where 


X(X)  = TO)  f d^A(x-x')V(x')A(x') 


(2.17) 


For  our  homogeneous  background  case,  where  the  Green's  function  A is 
given  by  Eq.  (2.5),  we  have 


X(x)  = — f d3x'  ^ eik(x' 

n J x'  I X-x'  I 


+|x-x  l-x) 


V(x') 


(2.18) 


In  concluding  this  section,  and  as  an  aside,  let  us  comment  on  the 
geometrical  optics  limit  of  this  expression.  This  limit  results  from  an 
evaluation  of  X(x)  by  the  method  of  stationary  phase.  Provided  that  the 
Fresnel  condition 


x < k L 


All  of  these  statements,  we  remind  the  reader,  are  subject  to  the  range 
constraints  under  which  the  supereikonal  expression  was  derived  in  the 
first  place. 


- •.•M&stigeprr*'' >‘i^<t‘" : .‘Ti1 


is  met,  the  stationary  phase  path  in  Eq.  (2.18)  is  the  straight  line  join- 
ing 0 to  x,  and  the  stationary  phase  value  of  X(x)  is  just 


X(x , 0 , 0)  = 


i r 

2k  / 

J n 


dx' V(x' ,0,0) 


(2.19) 


which  is  immediately  recognized  as  the  correct  geometrical  optics  expres 
sion  for  the  phase.  The  analogous  expression  for  the  amplitude  in 
geometrical  optics  is  obtained  by  keeping  the  second-order  transverse 
derivatives  in  V as  well. 
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III  FLUCTUATIONS  IN  A HOMOGENEOUS  ISOTROPIC  OCEAN 


I 


The  quantities  that  it  will  be  of  interest  to  compute  are  the 
statistical  averages  <X2(x)>  and  <|x(x)[2>.  (Wo  note  that  <X(x)>  = 0, 
of  course.)  These  are  connected  to  phase  and  amplitude  fluctuations 

2 99 

($  ) and  (|log  p/Pq|  ) = (A  ) by  the  relations 


(<l> 


(O  = 1/2  ( < I Xj  ) - Re<X“> 


) 


(3.1) 


4 

I 

fi 

1 

ift 

! 


and 


("■ 


) 


<A  > = 1/2  ( < | X|  > 4-  Rc<X  ) 


9 

We  will,  in  addition,  obtain  the  cross  correlation  Im  (x  ), 


(3.2) 


Insofar  as  the  sound-speed  fluctuations,  and  hence  the  fluctuations 
in  X,  are  gaussian,  the  pressure  fluctuations  arc  related  to  the  statisti- 
cal averages  as  well.  We  have 


/ \ 1/2 <X  > 

SP)  = PQe 


$ 

s 


, 2X  2 2<x“) 

\P  > = PQ  G 


(|p|  ) = !p0!  e 


21/2 <(x+X*)2> 


§ 

i 


where  we  write 


po(x) 


5 A (x)  = 


ikx 


4jtx 


I 


if 


* 

i 

1 

» 

m ■ 


tfi»i 


- KMttS  *■  - ■ ■ • . -,j 


as  the  received  pressure  from  the  unit  point  source  in  the  absence  of 


fluctuations. 


Eventually,  we  will  also  be  interested  in  correlations;  these  will 
involve  averages  such  as  (x&^XO^)  > , etc.,  but  we  shall  ignore  these 


for  now. 


Let  us  first  evaluate  <|x(x)|  >.  For  convenience,  we  shall  choose 
x to  lie  along  the  x axis,  so  that  x = (x,0,0).  From  the  definition, 

Eq . (2.18),  we  have 


<|x(x)|2>  = 


2 | | | -> 

y,  !x  - y,  |yJx  - y„ 


ikCy^lx-yJ-x)  -ik(y2+|x-y2|-x) 


:(l  ■ >2) 


(3.3) 


where  we  have  introduced  the  correlation  functi( 


“> 


C(YL  - y2)  s <V(y1)V(y2)> 


We  assume  C to  be  independent  of  (y^  +>  y^/2  for  this  case  of  a homo- 
geneous background. 


It  is  convenient  in  Eq . (3.3)  to  shift  to  relative  and  center-of- 


mass  coordinates.  We  define 


- y,  + y2 

y = yi  “ y2  ’ Y = — 


• S'! v ’■  * ■“ ' ' *0  » " • 

•Tf  % J.KT’1  ?,r 
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Then,  if  we  assume  that  C(y)  cuts  off  lor  values  of  y ^ L where  L « x, 
we  may  expand  in  y/Y.  Thus,  Eq.  (3.3)  becomes 


<|x(x)|2>  = 


■ (ij  {■ 


2 _ 2 
Y Y - x 


/■' 


y C(y) 


exp  ik  y • - (x  - Y)] 


(3.4) 


Here  Y and  x - Y stand  for  unit  vectors  in  the  direction  of  Y and  x - Y, 
respectively,  and  we  have  written  |y  ± y/2 1 » Y,  |x  - Y ± y/2 1 |x  - y| 
in  the  geometrical  factors  multiplying  the  exponentials.  This  approxi- 
mation introduces  a negligible  error. 


3-> 


The  integral  over  d Y may  now  be  evaluated  by  stationary  phase.  The 


stationary  phase  path  is  the  straight  line  joining  0 to  x,  and  the  result 
is 


/ 


< | X (x)  | 2)  = / dyn  C(yM  ,0) 

4k 


(3.5) 


whore  y^  refers  to  tile  component  of  y in  the  direction  parallel  to  x. 
Introducing  the  Fourier  transform  of  the  correlation  function 


/->  -> 

3->  -iq.r  -> 
d r e C(r) 


(3.6) 


pei'mits  us  to  rewrite  Eq . (3.5)  in  the  sometimes  more  convenient  for 


•m 


<l«sl2>  -fe) 


2->  -> 

q J.  c(0,q±) 


(3.7) 


whore  "x"  refers  to  the  directions  perpendicular  to  x. 
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At  this  point  it  is  convenient  to  express  C (y^ , y L)  in  terms  of  its 
iourier  transform,  as  given  by  Eq.  (3.6).  The  integral  over  cly ^ d2yx 
can  then  be  carried  out,  and  we  finally  obtain  the  relatively  simple 
expression 


(x (x)  / = - [ — 
\ 4nk 


2q±  C(0,q  ) 


l 


R±  (s-x)s 


ds  e 


(3.11) 


Equations  (3.7)  and  3.11)  constitute  our  central  results.  They  ex- 
press the  quantities  of  interest  as  integrals  along  unperturbed  ray  paths 
(in  this  case  straight  linos)  of  the  Fourier  transform  of  the  correlation 
function  C(q)  times  rather  simple  geometrical  factors.  As  we  shall  see 
later,  entirely  parallel  expressions  obtain  in  the  more  difficult  case  of 
an  inhomogeneous  background  medium. 

The  expression  for  <|X(x)|  ),  Eq.  (3.7),  is  precisely  the  same  result 
for  this  quantity  obtained  by  using  geometrical  optics  to  compute  X(x)  it- 
self, and  then  calculating  (|x(x)|  ) from  this,  [This  is  easily  seen  by 
referring  back  to  Eq.  (2.19).]  In  contrast,  Eq.  (3.11)  is  not  what  one 
obtains  for  <X(x)2)  from  geometrical  optics.  Geometrical  optics  for  this 
quantity  is  recovered  if  one  expands  the  exponential  in  Eq,  (3,11),  a 
procedure  that  evidently  is  valid  only  if 


Ql  (s  - x)s 


« 1 


Since  qx  ~ 1/L  and  s,  x - s ~ x,  this  condition  can  be  more  familiarly 
written  as 


1 

■A 


& 

&' 


.n 

1 


..i 


m 

ym 


M 

$ 


13 


x « k L2 

which  we  recognise  us  the  Fresnel  condition  under  which  the  geometrical 
optics  approximation  lor  X(x)  itself  was  valid  in  the  first  place. 

lhus,  Eq . (3.11)  constitutes  an  improvement  over  geometrical  optics, 
while  Eq.  (3.7)  coincides  with  geometrical  optics.  Conversely,  geometri- 
cal optics  for  (|x|  ) is  valid  out  to  a very  large  range,  while  geometri- 
cal optics  for  (x2)  is  valid  only  within  the  range  x < k L2. 

It  is  of  interest  to  study  Eq.  (3.11)  in  the  limit  of  very  long 
range.  As  x - ® , the  integral  over  ds  can  bo  approximately  evaluated, 
and  we  find1 


<X2(x)> 


iC  (0) 
8rck 


(V  + log  4kx  - iit/2) 


(3.12) 


where  V — 0. 577,., is  Euler's  constant.  Thus,  for  small  x satisfying 
the  Fresnel  condition,  we  have 


<X2(x)> 


i 


6(0,^) 


[X  IT  X + --J  • (3.13) 

and  for  very  large  x we  have  (x2(x))  rvJ  i log  x as  given  by  Eq,  (3,12). 
Between  these  two  extremes,  of  geometrical  optics  and  of  very  long  ranges, 
Eq.  (3,11)  provides  a smooth  interpolation. 

p 

Equation  (3.7),  on  the  other  hand,  gives  (|x|  ) for  all  values  of  x, 
largo  and  small,  and  simply  says  that  (|x|2)  is  proportional  to  the  range 
x everywhere. 


14 
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IV  RYTOV'S  METHOD  IN  AN  INHOMOGENEOUS  OCEAN 


Now  let  us  turn  to  the  effects  of  the  sound  channel.  That  is,  we 

must  replace  the  nonfluctuating  sound  speed  c in  the  homogeneous  case  by 

a (specified)  function  of  position  c(x).  In  fact,  for  the  ocean, 

c(x)  = c ( z ) is  a function  of  depth  only.  A reasonably  explicit  form  for 

3 

c(z)  is  the  relatively  simple  expression 


c(z)  = c [1  + e (e“  + T)  - l)] 

A 


where  T|  = (z  - z )/£  B,  where  z is  the  sound-channel  depth,  and  e and  B 
A A 

aie  paiameters.  We  shall,  however,  write  our  formulas  for  a general  c(x) 
until  the  time  comes  to  make  explicit  numerical  estimates. 

The  wave  equation  for  the  pressure,  which  is  our  starting  point,  now 
becomes  altered  from  Eq.  (2.1)  to  the  equation 


[V  + k^(x)]p(x)  = V(x)p(x) 


(4.1) 


still  with  the  same  boundary  condition  that 


as  x " 0,  corresponding  to  an  isotropic  point  source  at  the  origin,  where 
now  k(x)  = in/c(x) . 

We  must  first  study  the  nonfluctuating  part  of  the  problem,  to  evalu- 
ate the  Green's  function  in  the  presence  of  the  sound  channel.  This 
satisfies 


aaswwasmaaKaawr^  i'-W*- 


[V2  + ^(xJjACx^y)  = 63(x-y)  . (4.2) 

Note  that  it  is  no  longer  a function  only  of  x-y  as  it  was  in  the  homo- 
geneous case.  We  shall  assume  that  geometrical  optics  provides  a good 

approximation  to  the  nonfluctuating  sound-channel  problem.  This  moans 

— > , 
that  we  can  represent  A(x,y)  as  a sum  of  contributions  from  each  ray  join- 

— > “> 

ing  x and  y.  To  be  specific,  we  may  write 


n(x,y) 

A(x,y)  - ^ A.(x,y")  (4.3) 

i=l 


where  n(x,y)  is  the  number  of  rays  and  A is  the  contribution  of  the  i 
ray.  We  have,  in  particular  for  rays  joining  the  origin  and  x, 


Aa (x,0) 

where  ds  is  an  element  of  path  length  along  the  ray,  x. (s)  is  the  i ray 
-> 

joining  0 to  x,  and  K.  is  a normalization  factor. 

Now  when  the  fluctuations  are  turned  on,  the  signals  traveling  on 
each  of  the  rays  joining  the  origin  to  the  point  of  observation  x are 
subject  to  small-angle  scatterings  by  the  perturbing  potential  V(x).  The 
signals  are  thus  deflected  slightly  from  the  undisturbed  rays  by  each  in- 
tern- tion  with  V.  The  repeated  action  of  V thus  produces,  on  each  ray,  a 
sort  of  random  walk  of  the  signal  away  from  the  original  ray.  When  we 
average  over  an  ensemble  of  perturbations  V,  the  disturbed  signals  will 
fill  up  a tube  surrounding  the  undisturbed  ray.  Provided  that  these  tubes 


—) 

K (x,0)  exp 

l 


ds  k(xi(s) ) 


i = 1. . ,n(x) , 


(4.4) 
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around  each  of  the  original  rays  do  not-  overlap,  the  rccicved  pressure 
will  be  a sum  of  contributions  from  each  ray  tube. 

We  may  estimate  the  radius  of  a ray  tube  as  follows.  The  mean  free 
path  d between  interactions  of  the  signal  traveling  along  a given  ray  with 
the  perturbing  potential  V is  of  the  order  of  kc/6c.  Hence,  over  a range 
x the  number  of  scatterings  is  n = x/d.  The  average  deflection  angle  due 
to  each  scattering  is  of  the  order  of  1/kL^  vertically  and  1/kL^  horizon- 
tally. where  L and  L are  the  vertical  and  horizontal  correlation  lengths 
of  the  sound-speed  fluctuations.  Since  the  process  is  a random  walk,  the 
net  displacement  due  to  n collisions  is  proportional  to  n,  and  hence  the 
vertical  and  horizontal  extents  of  the  tube  are,  roughly, 


j » » w 


Lot  us  assume  that  the  vertical  extent  of  the  tubes  is  small  enough 
so  that  the  tubes  remain  distinct.  Then  the  pressure  at  x is  the  sum  of 
contributions  from  each  tube, 


P . (x) 
1 


(4.5) 


where  n(x)  is  the  number  of  unperturbed  rays  joining  the  source  to  the 
point  x.  We  shall  bo  interested  in  (x) , 


■ 1'  " . -'s  ifcV*  \ \ . 

*i»  1 1 ’*  ft  JV  f- 

-ft.  * ji?  'A1*-'"*;;.  x 


fJV 


Wc  n°te  that  pi(x)  is  the  Pressure  that  would  be  received  at  x if 
the  source  were  not  isotropic,  but  rather  emitted  all  its  energy  in  the 
direction  of  the  iUl  unperturbed  ray.  Thus  p. (x)  must  satisfy  the  wave 
equation  (4.1)  but  with  an  anisotropic  boundary  condition  that  itself 
depends  on  x.  To  make  this  more  precise,  let  us  deirne  p.(y;x)  to  be 
the  pressure  at  y from  a source  at  the  origin  that  emits  only  within  a 
small  solid  angle*  around  the  direction  of  the  ith  unperturbed  ray  join- 
ing the  origin  to  3.  Thus  p. (3)  = P.(x;x),  and  furthermore  p.(y,x) 
vanishes  unless  y is  inside  the  ith  ray  tube.  Then 


[V-> 


y + k (3' >3  p±  Cy  ;x>  = v(y)p.(y;x),  i = l...n(x) 


r*  -» 


(4.6) 


Tn  analogy  with  this  definition  of  p^yjx),  we  may  also  define  an 
unperturbed"  Green's  function  A.(y;x,0),  i = l...n(x),  to  satisfy 


tVy  + k2(y)]Ai (y;x,0)  = 0 


(4.7) 


again  with  the  same  boundary  condition.  This  function,  also,  will  vanish 
except  when  y is  near  the  i th  unperturbed  ray. 


We  may  now  directly  derive  the  analogue  of  Eq.  (2.15)  by  computing 
the  quantity  log  p. (y ;x)/A . (y ;x,o)  in  perturbation  theory,  and  using 
Eq.  (4.2).  We  find,  setting  y = x,  that 


P.(x)  = A.(x;x,0)e 


~rt  f 

_>  A . (x;x,0)  /th 

^ n\„  1 r: 


3 - >,  -*  -»  -»  -* 

d x A(x,x  )V(x'  )A(x'  ,0) 

ay  tube 


(4.9) 


We  have  here  replaced  A.(x;x,0)  simply  by  A.(x,0),  Equation  (4.9)  is 
evidently  the  generalization  of  the  Rytov  formula  (2.15)  to  the  situation 
of  an  inhomogeneous  background  and  many  rays.  The  expression  clearly 
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fails  if  the  range  is  so  large  that  thr  ray  tubes  overlap;  otherwise  the 
validity  conditions  are  the  same  as  those  in  the  homogcncous-baekgi’ound 
case . 
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V FLUCTUATIONS  IN  AN  INHOMOGENEOUS  OCEAN 


In  this  section  we  shall  use  Eq.  (4.9)  to  calculate  the  various 
averages  of  interest  for  the  contribution  of  a single  ray  tube  to  the 
pressure  in  the  presence  of  the  sound  channel.  We  shall,  for  simplicity, 
drop  the  index  i,  though  we  should  keep  in  mind  that  when  there  are 
several  unperturbed  raypaths  their  contributions  are  to  be  added  to  obtain 
the  total  pressure.  Our  interest,  then  will  be  in  the  statistical  fluc- 
tuations of  the  contributions  of  a single  ray,  or  rather  a single  ray  tube. 

As  in  the  homogeneous  ease,  we  define 


X(x) 


-±-fj 

i ( x , 0 ) J 


y A(x,y)V(y)A(y,0) 


(5.1) 


2 i 1 2 

and  we  wish  to  compute  <X  > and  <|x|  ).  We  recall,  from  Seetion  IV,  that 
assuming  geometrical  optics  to  be  a valid  approximation  for  the  nonfluc- 
tuating background  permits  us  to  write  the  Green's  function  as 


a <s,?>  = Ka,WkS(i 


(5.2) 


where 


-»  -♦  r y 

kS(x,y)  = / ds  k(x(s) ) 


(5.3) 


and  where  the  normalization  faetor  is 
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Flore,  x refers  to  directions  perpendicular  to  the  ray.  (An  excellent 
approximation  to  tliis,  for  our  purposes,  is  to  write  simply  K(x,y)  = 

( 4 Tt  I x— y | ) , as  in  the  liomogeneous-ocean  case;  we  need  to  Foe  careful 
about  deviations  from  homogeneity  only  in  the  phases.)  In  Eq.  (5.3) 
the  line  integral  is  along  the  ray  of  interest  joining  the  points  x 


and  y. 


To  repeat  our  earlier  calculations  requires  us  to  introduce  the 
correlation  function  (v(y^) V(y^))  . In  the  homogeneous  case,  this  quan- 
tity depended  only  on  the  separation  y - y Now,  however,  because  of 

tlie  background  inhomogeneities,  it  will  also  depend  on  (y  + y ) /2 

X 2 

(actually  it  will  depend  only  on  the  mean  depth  (z  + z )/2  since  the 

X 2 

inhomogeneities  depend  only  on  depth).  Thus  we  must  now  define  the 
correlation  function  by 


V<yi’y2)  * C Pi  " V 


y + y 
*1  2 


As  before,  let  us  look  first  at  <|x|2>.  We  have 


<lx(x)i2>  = o)  c(-t) 

J K(x, 0)  J 


exp  i k [ S ( x ,Y  + y/2)  - S(x,Y  - y/2) 


+ S(Y  + y/2 ,0)  - S(Y  - y/2,0)] 


(5.4) 


and  we  must  keep  in  mind  that  we  are  to  integrate  only  over  the  ray  tube 
suiiounding  the  unpci turbed  ray  of  interest.  In  the  homogeneous-background 
case  we  expanded  the  exponent:  in  powers  of  y,  because  C(y)  vanished  for 
large  |y|.  We  may  do  the  same  here.  Thus, 


( 


x(x>  1 2>  = f ak^v)k(v,°)  rci3_.  cGt 


Y) 


K(x,0) 


exp  i k y • V [S(x,Y)  + S(Y,0)] 


(5.5) 


The  integral  on  d Y is  again  to  be  evaluated  by  stationary  phase.  The 
stationary  phase  path  is  evidently  the  unperturbed  ray  joining  the  origin 
to  the  observation  point  x.  Hence,  we  may  write 


<|x(x)|2> 


2_, 

d qA(s)  C(qx (s) F Y(s))  , 


(5.6) 


in  complete  parallel  to  the  homogeneous  case.  Here 
ds  is  along  the  unperturbed  ray,  qx(s)  refers  to  the 
pendicular  to  the  ray  at  s,  Y(s)  is  a point  on  the  r 

C(q,Y)  s / d3y  eiC1*yc(y,Y) 


the  line  integral  on 
component  of  cf  per- 
ay  at  s,  and 


(5.7) 


Next  we  turn  to 


<X2>: 
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. -,  2.  f 3-.  K(x , Y)  K(Y , 0)  f 3-*  -.  - 

<X(x)  ) = / d Y — / d y C(y ,Y) 

J K(x ,0)  J 


exp  i k 


[s(x,Y  + 


y/2)  + S(x,Y  - y/2) 


+ S (Y  + y/2 ,0)  + S(Y  - y/2,0) 


- 2 S(x 


,0)] 


(5.8) 


Now  when  we  expand  the  exponent  in  powers  of  y the  linear  terms  vanish,  so 
that  we  have 


. _ 2.  C 3-  K(x, Y)K(Y,0)  2ik (S (x , Y)+S (Y , 0) -S (x , 0) ) 

(X(x)  ) = I d Y e 

J K(x,0) 


J d3y  C(y , Y)  exp  ik/4 


(5.9) 


where  we  define 


A.  .(Y) 
ij 


J_  _L  [sU.V)  ♦ SCY,0)] 
i J 


(5.10) 


Evaluation  of  the  integral  on  d3Y  by  stationary  phase  again  selects  as  the 

— > 

stationary  phase  path  the  unperturbed  ray  joining  0 to  x.  The  integral  on 
d y can  then  be  done  by  introducing  the  Fourier  transform  C(q,Y)  as  in 
Eq,  (5,7).  Finally,  we  obtain 


<X(x)  > = 


-)/ 


q±(s) 


C(q  (s),  Y(s))e 


r q,  . (s)q.  . (s) A (Y(s) ) 
k 1 j.  1j  ij 


(5. 11) 


The  notation  is  as  in  Eq.  (5.6),  and  the  result  is  again  in  complete  analogy 
to  the  homogeneous  case. 


Most  of  the  comments  we  made  in  Section  III  concerning  the  results  in 
the  homogeneous  background  apply  here  as  well.  The  expression  for  (|x|  ) 

is  again  just  that  obtained  in  the  geometrical  optics  approximation,  but 

2 2 
that  for  (X  ) is  not.  Geometrical  optics  for  (X  ) is  valid  provided  that 


7 qL.  (s)  q.  .(s)  A 1 ( Y ( s ) ) . . « 1 

k xi  ±J  ij 


(5.12) 


This  is  the  analogue  of  the  Fresnel  condition. 

In  the  homogeneous  case,  this  condition  boiled  down  to  the  requirement 


x < k L 


In  the  presence  of  a sound  channel,  the  restriction  (5.12)  on  the  range  is 
less  severe;  the  condition  (5.12)  reduces  approximately  to 


x x tan  9 _ 

4-  « 1 

2 2 
k L k L 

H V 


(5.13) 


where  tan  9 is  the  maximum  ray  inclination  to  the  horizontal.  Since 
is  much  larger  than  L^,  and  since  tan  9 is  small,  this  restriction  is 


^ '.T  " ' V'ss’t  ‘ 
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easier  to  meet  than  x/k  L * « 1.  Thus,  in  the  presence  of  a sound 
channel,  geometrical  optics  should  be  valid  to  a greater  range  than 
would  be  the  case  with  a uniform  background  sound  speed. 


Wo  have  presented  relatively  simple  formulas  for  phase  and  amplitude 
fluctuations  of  pressure  signals  in  the  ocean  in  the  presence  of  a speci- 
fied spectrum  of  sound-speed  fluctuations  for  both  a homogeneous  ocean 
and  one  with  a sound  channel.  These  formulas  reduce  to  the  geometrical 
optics  approximation  for  short  ranges,  but  for  long  ranges  they  give  re- 
sults far  less  divergent  with  range  than  does  geometrical  optics. 


Foi  the  case  of  a sound  channel,  where  there  are  in  general  many  rays, 
our  results  apply  to  fluctuations  in  the  contribution  of  a single  ray  to 
the  received  pressure;  thus  the  approximation  is  limited  to  ranges  at  which 

■ 

lays  are  still  separable.  Fluctuations  in  the  total  received  pressure  are 
much  greater,  due  to  interference  between  different  rays,  and  are  insensi- 

* 

tive  to  details  of  the  fluctuations  in  a single  ray,4 

The  next  step  to  be  undertaken  is  to  use  the  results  outlined  here, 
together  with  a semi-empirical  fluctuation  spectrum,  to  make  numerical 
estimates  of  the  fluctuations  for  the  purpose  of  comparison  with  experi- 


ments. This  will  be  described  elsewhere. 

! 
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